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ABSTRACT 



Aims. The interacting dark matter (IDM) scenario allows for the acceleration of the Universe without Dark Energy. 

Methods. We constrain the IDM model by using the newly revised observational data including H(z) data and Union2 SNe la via the 

Markov chain Monte Carlo method. 

Results. When mimicking the ACDM model, we obtain a more stringent upper limit to the effective annihilation term at kC\ ~ 
10~ 3 ' 4 Gyr~', and a tighter lower limit to the relevant mass of Dark Matter particles at M x ~ 10~ 8 - 6 Gev. When mimicking the wCDM 
model, we find that the effective equation of state of IDM is consistent with the concordance ACDM model and appears to be most 
consistent with the effective phantom model with a constant EoS for which w < — 1. 
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1. Introduction 

Recent observations of type la supernovae (SNe la, Riess et al. 
1998; Perlmutter et al. 1999) have predicted that our present uni- 
verse is passing through an accelerated phase of expansion pro- 
ceeded by a period of deceleration. A new type of matter with 
negative pressure, which is popularly known as dark energy, has 
been proposed to explain the present phase of acceleration. The 
most simple dark energy candidate, the cosmological constant 
(ACDM model), though known to be consistent with various ob- 
servations such as SNe la, the galaxy cluster gas mass fraction 
data (Wilson et al. 2006; Davis et al. 2007; Allen et al. 2008), 
and the CMB temperature and polarization anisotropics (Jassal 
et al. 2010), is always affected by the coincidence problem. Until 
now, many other dark energy models have been brought for- 
ward to explain this comic acceleration such as the scalar fields 
with a dynamical equation of state [e.g., quintessence (Peebles 
& Ratra 1988a, 1988b; Caldwell, Dave, & Steinhardt 1998), 
phantom (Caldwell 2002), k-essence (Armendariz-Picon et al. 

2001) , quintom (Feng et al. 2005; Guo et al. 2005; Liang et al. 
2009)], the Chaplygin gas (Kamenshchik et al. 2001; Bento et al. 

2002) , holographic dark energy (Cohen 1999; Li 2004), and so 
on. Many alternatives to dark energy in which gravity is modified 
have been proposed as a possible explanation of the acceleration 
[e.g., the braneworld models (Dvali, Gabadadze, & Porrati 2000; 
Zhu & Alcaniz 2005), the Cardassian expansion model (Freese 
& Lewis 2002, Zhu & Fujimoto 2002)], as well as the f(R) the- 
ory (Capozziello & Fang 2002; Carroll et al. 2004) and the f(T) 
theory (Bengochea & Ferraro 2009; Wu & Yu 2010). 

It has been shown that the dark matter self-interactions could 
provide the accelerated expansion of the Universe without any 
dark energy component (Zimdahl et al. 2001; Balakin et al. 

2003) . In the framework of the Boltzmann formalism, if there 
is a disequilibrium between the dark matter particle creation 
and annihilation processes, an effective source term with neg- 
ative pressure could be created. Basilakos & Plionis (2009) in- 



vestigated the circumstances under which the analytical solution 
space within the framework of the interacting dark matter (IDM) 
scenario allows for a late accelerated phase of the Universe, and 
find that the effective annihilation term of the simplest IDM 
model is quite small by using the nine observational H(z) data 
points of Simon et al. (2005). The gravitational matter creation 
model that is fully dominated by cold dark matter (Lima et al. 
2008; Basilakos & Lima 2010) is mathematically equivalent to 
one case of the IDM models (Basilakos & Plionis 2009). 

In this paper, we use the newly revised H(z) data (Stern et 
al. 2010; Gaztanaga et al. 2009) and the Union2 set of 557 SNe 
la (Amanullah et al. 2010) to constrain the relevant IDM models 
(Basilakos & Plionis 2009) by using the Markov chain Monte 
Carlo (MCMC) method. This paper is organized as follows. In 
Sect. |2l we briefly indicate the basic equations of the IDM mod- 
els. Observational data including H(z) and SNe la are given in 
Sect. [3] In Sect.|U MCMC constraint results from different com- 
bined data sets are illustrated. Finally, we summarize our main 
conclusions in Sect. [5] 

2. The basic equations of the IDM models 

We assume that the IDM density obeys the collisional 
Boltzmann equation (Kolb & Turner 1990) 

(1) 



p + 3Hp + k P 2 - ¥ = 



where *F is the rate of creation of DM particle pairs and the an- 
nihilation parameter k = (cru)/M x (where M x is the mass of the 
DM particle, cr is the cross-section for annihilation, and u is the 
mean particle velocity). Compared to the usual fluid equation, 
the effective pressure term is 



Peff = {k P 2 - V)/3H . 



(2) 



When the IDM particle creation term is larger than the annihi- 
lation term (*p 2 - W < 0), IDM may serve as a negative pres- 
sure source in the global dynamics of the Universe (Zimdahl et 
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al. 2001; Balakin et al. 2003). Basilakos & Plionis (2009) phe- 
nomenologically identified two functional forms for which the 
previous Boltzmann equation can be solved analytically, only 
one of which is of interest since it indicates the dependence 
of the scale factor on a "oc a -3 ". We refer to appendix B in 
Basilakos & Plionis (2009) for details. We assume that 



y(a) = aH(a)R(a) ■■ 



C l (n + 3)a"H(a) + KC 2 l a 2n 



And the total energy density is 

f \ n n , -3 F ( a ) 
p(a) — L\a + a 



[C 2 - f 1 a x- 3 f(x)F(x)dx] 



(3) 



(4) 



where n, C\, and C2 are the corresponding constants of the 
problem (kC\ in the unit of Gyr~'), and the kernel function 
F(a) = exp[-2/cCi x"' 1 /H(x)dx\. The first term on the right 
side of Eq.© obviously corresponds to the residual matter cre- 
ation that results from the possible disequilibrium between the 
IDM particle creation and annihilation processes, while the sec- 
ond term can be viewed as the energy density of the self IDM 
particles that are dominated by the annihilation process. 

2. 1 . Model 1 : Mimicking the KCDM Model 

If n = 0, the global density evolution can be transformed as 



p(a) — C\ + a 



e -2 K Ci(t-ta) 

[c 2 -«z(0] ' 



(5) 



where Z(f) = f flrV 2 ^ 1 ( '-' o) (Basilakos & Plionis 2009). At 
the present epoch, the density evolves according to p(a) C\ + 
a~ 3 /C2, which is approximately equivalent to the corresponding 
evolution in the ACDM model, in which the C\ term resembles 
the cosmological constant term (pa) and the I/C2 term resem- 
bles the density of matter (p,„). The Hubble parameter can be 
written as (Basilakos & Plionis 2009) 



— I =Qi j0 + C22,oa 



e 



-2Kd(t-t ) 



l + KCi(fi 2 ,o/Oi,o)Z(0 



(6) 



where Qi,o = SnGCi /3H 2 and Q 2 ,o = 87rG/3// 2 C 2 , which relate 
to £2a and Q m o in the ACDM model, respectively. The mass of 
the DM particle can also be related to the range of kC\ (Basilakos 
& Plionis 2009) 



M x 



1.205 x 10~ 12 (cru) 



kC] 



10 



-22 



GeV 



(7) 



2.2. Model 2: Mimicking the wCDM Model 

When the annihilation term is negligible (k = 0) and the parti- 
cle creation term dominates, it is straightforward to obtain the 
evolution of the total energy density (Basilakos & Plionis 2009) 



p(a) = Cifl" +(1/C 2 K 



(8) 



The conditions n > -2 implies that the IDM fluid has an inflec- 
tion point, and the current model acts as the wCDM model with 
a constant EoS with widm = -1 - n/3, but with a totally differ- 
ent intrinsic nature. This situation is mathematically equivalent 
to the gravitational DM particle creation process within the con- 
text of non-equilibrium thermodynamics (Lima et al. 2008). The 
Hubble parameter is now given by 



-3 



(9) 



where Q 2 , = 8ttG(C 2 - Ci)/3// 2 and S2 1>0 = %nGCi/3H 2 , re- 
spectively (Basilakos & Plionis 2009). 



3. Observational data 

To constrain the relevant IDM models (Basilakos & Plionis 
2009), we use the newly revised H(z) data (Stern et al. 2010; 
Gaztanaga et al. 2009) and the Union2 set of 557 SNe la 
(Amanullah et al. 2010). 



3.1. The observational H(z) data 

It is known that the Hubble parameter H(z) depends on the dif- 
ferential age as a function of redshift z in the form 



H(z) 



1 dz 
1 +z dt' 



(10) 



which provides a direct measurement on H(z) based on dz/dt. 
Jimenez et al. (2003) demonstrated the feasibility of this method 
by applying it to a z ~ sample. Moreover, compared with other 
observational data, it is more rewarding to investigate the obser- 
vational H(z) data directly, because it can take the fine structure 
of H{z) into consideration and then use the important informa- 
tion that this structure provides. 

By using the differential ages of passively evolving galax- 
ies determined from the Gemini Deep Deep Survey (GDDS) 
(Abraham et al. 2004) and archival data (Treu et al. 2001, 
2002; Nolan et al. 2003a, 2003b), Simon et al. (2005) deter- 
mined nine values of the Hubble parameter H{z) in the range 
< z < 1.8, which have been used to constrain the param- 
eters of various dark energy models (Samushia & Ratra 2006; 
Wei & Zhang 2007; Wu & Yu 2007a, 2007b; Zhang & Zhu 
2007; Kurek & Szydlowski 2007; Lazkoz & Majerotto 2007; 
Sen & Schen-er 2008; Yi & Zhang 2007; Wan et al. 2007; 
Xu et al. 2008; Zhai, Wan, & Zhang 2010). The H(z) data 
at 1 1 different redshifts were determined from the differential 
ages of red-envelope galaxies (Stern et al. 2010), and two more 
Hubble parameter data points H(z = 0.24) = 79.69 ± 4.61 and 
H(z = 0.43) = 86.45 ± 5.96 were obtained by Gaztanaga et 
al. (2009) from observations of BAO peaks [for a review of the 
observational H(z) data, see Zhang & Ma (2010)]. Studies us- 
ing these newly H(z) data for cosmological constraint include 
Gong et al. (2010), Liang et al. (2010a), Liang & Zhu (2010), 
Cao & Liang (2010), Ma & Zhang (2011), and Xu & Wang 
(2010c). We emphasize the use of the Hubble constant Hq in 
our analysis. Many previous have determined its present value. 
Freedman et al. (2001) presented the final results of the Hubble 
Space Telescope (HST) key project that measured the Hubble 
constant Hq -12 + 8kms _1 Mpc _1 , Gott et al. (2001) and Chen 
et al. (2003) proposed that Ho = 68kms~ 1 Mpc~ 1 was a more 
likely value, and Tammann et al.(2008) obtained Hq = 62.3 + 
1.3kms Mpc 1 from 28 independently calibrated Cepheids and 
the distant, Cepheid-calibrated SNe la. More recently, Riess et 
al. (2009) determined Hq = 74.2 ± 3.6kms _I Mpc _1 by com- 
bining the observations of 240 Galactic Cepheid variables using 
HST. In this work, we follow Basilakos & Plionis (2009) and 
adopt Hq = 72 + 8kms _1 Mpc _1 . The observational H(z) data are 
given in Table [TJ In this case, x 2 can be defined as 



14 



-z 



(H - H b s ) 



(11) 
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0.1 


0.17 


0.24 


0.27 


0.4 


0.43 


0.48 


0.88 


0.9 


1.3 


1.43 


1.53 


1.75 


H(z) 


72 


69 


83 


79.69 


77 


95 


86.45 


97 


90 


117 


168 


177 


140 


202 


lcr uncertainty 


±8 


±12 


+8 


+4.61 


±14 


±17 


±5.96 


±60 


±40 


±23 


±17 


±18 


±14 


±40 



Table 1. The observational H(z) data (Stern et al. 2010; Gaztanaga et al. 2009; Freedman et al. 2001) 



3.2. The observational SNe la data 

SNe la have long been used as "standard candles". It is com- 
monly believed that measuring both their redshifts and apparent 
peak fluxes gives a direct measurement of their luminosity dis- 
tances, thus SNe la provide the strongest constraints on the cos- 
mological parameters (Riess et al. 2004, 2007; Astier et al. 2006; 
Davis et al. 2007; Wood-Vasey et al. 2007; Kowalski et al. 2008; 
Hicken et al. 2009; Chen et al. 2010). The present analysis uses 
the Union2 (557 sample) data set of the Supernova Cosmology 
project covering a redshift range 0.015 < z < 1.4 (Amanullah et 
al. 2010), which was used to constrain cosmological models in 
Wei (2010), Liang et al. (2010a, 2010b), Liang et al. (201 1), Xu 
& Wang (2010a, 2010b), and Wang et al. (2010). 

In the calculation of the likelihood from SNe la, we 
marginalize the nuisance parameter by minimizing (Di Pietro & 
Claeskens 2003), 

^^-fWl)' (12) 

where A = (// data - /J. th ) 2 /cr 2 , B = £/ 5? (^ data - j"" 1 )/^- , 
C = Jip 7 l/of> an d the distance modulus is /j = 
5 log^/Mpc) + 25, with the lcr uncertainty cr, from the ob- 
servations of SNe la and the luminosity distance dh as a function 
of redshift z, 

4. Constraint results 

We define the total likelihood to be the product of the separate 
likelihoods of the two cosmological probes, in other words, 

*Ld=*ff +*SNe- (I 4 ) 

To determine the best-fit cosmological parameters, we use a x 2 
minimization and the 68.3% and 95.4% confidence levels are 
defined by the conventional two-parameters x 2 levels 2.30 and 
6.17, respectively. We perform a global fitting to determine the 
cosmological parameters using the Markov chain Monte Carlo 
(MCMC) method and our MCMC code is based on the publicly 
available package COSMOMC (Lewis & Bridle 2002). 

4.1. Model 1: Mimicking the ACDM Model 

In this model, there are two free parameters: Q2,o and kC\ (or 
M x ) for a flat background (Q^o + Qift = 1)- In this paper, 
we follow Basilakos & Plionis (2009) and adopt fo = l/#o = 
13.6 Gyrs (roughly the age of the universe in the corresponding 
ACDM cosmology). The constraint results from different data 
combinations are shown in Fig.Q]-[3]and summarized in Table|2] 
In Fig.Q] we show the results from the H(z) data. The best- 
fit solution corresponds to fit Q2,o = 0.270^}^ and log(A:Ci • 
Gyr) w -6.85 with an upper limit log(/cCi ■ Gyr) » -3. This 
effective annihilation term is obviously still unconstrained to- 
wards lower values. Correspondingly, the best-fit value of M x 



Model 1 


H(z) 


SNe 


tf(z)+SNe 


-3 

X min 
£^2,0 

log(*Ci ■ Gyr) 
logM v /Gev 


9.86 

270 +0 - 044 

"•■^ '"-0.044 

-6.85(< -3) 
-5.15(> -9) 


544.13 

q 271+0 033 
KJ - CI 1 -0.031 

-4.42(< -3.5) 
-7.58(> -8.5) 


554.34 

f) 777+0.028 
"■^ -0.029 

-5.35(< -3.4) 
-6.65(> -8.6) 



Table 2. Summarizing the results of parameter constraints from 
model 1 considered in this work. We adopt fo = 1 1 Ho = 
13.6Gyrs in the computations. 




Fig. 1. The likelihood contours at the 68.3% and 95.4% confi- 
dence levels in the £^2,0 _ kC\ and Q2,o - M x planes provided by 
fitting model 1 to the H(z) data. 




Fig. 2. The likelihood contours at the 68.3% and 95.4% confi- 
dence levels in the £~22,o _ kC\ and £22,o - M x planes provided by 
fitting model 1 to the SNe la data. 



is logM A /Gev w -5.15 with a relatively stringent lower limit 
of logM A /Gev « -9. These results are consistent with previous 
work that uses fewer observational H(z) data 02,0 = 0.3^ ^and 
log(^d ■ Gyr) * -9.3 (Basilakos & Plionis 2009). In Fig.Sl we 
show the constraint results from SNe la with the best-fit model 
parameters Q 2 , = 0.27 l+J}^ and log(*Ci • Gyr) a -4.42. A 
more stringent upper limit is obtained at around log(/cCi • Gyr) » 
-3.5. The best-fit DM particle mass is logM x /Gev ~ -7.58 
with a relatively stringent lower limit logM^/Gev » -8.5. To 
obtain a tighter constraint on the model parameters, we com- 
bine the H(z) and SNe la data, and the results are shown in 
Fig. [3] The best-fit model parameters are £22,0 = 0.272^q q|| and 
log(A:Ci ■ Gyr) =s -5.35 with a much more stringent upper limit 
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Fig. 3. The likelihood contours at the 68.3% and 95.4% confi- 
dence levels in the ^2,0 _ kC\ and Q2 ; o - M x planes provided by 
fitting model 1 to the k(z)+SNe data. 




Fig. 4. The likelihood contours at the 68.3% and 95.4% confi- 
dence levels in the Q2,o - n and ^2,0 - widm pl anes provided by 
fitting model 2 to the H(z) data. 

log(/cCi ■ Gyr) * -3.4. Moreover, the best-fit value of M x is 
logM x /Gev w -6.65 with a relatively stringent lower limit of 
logM t /Gev « -8.6 at 2<x. Since M x is unbound at small val- 
ues, it is consistent with currently accepted lower bounds of M x 
(lOGeV) (see Cirelli et al.(2009) and references therein). 

4.2. Model 2: Mimicking the wCDM model 

In this case, there are two free parameters: Q.2,0 and n (or 
widm)- The constraint results from different data combinations 
are shown in Fig. [4]- Fig. [6] and summarized in Table [3] 



Model 2 


H(z) 


SNe 


H(z)+SNe 


X 1 

mm 
^2,0 

n 

WlDM 


9.40 

fl 977+0 098 
K '- z -' -0.097 

70 +1 54 

'"-1.54 
1 9^+0.51 


544.12 
0.288^ 

-1.06!$ 


554.39 
0.28« 
0.22:g 
-1.07+^ 



Table 3. Summarizing the results of constraint on parameters 
from model 2. 



For the H(z) data and the results shown in Fig. [4] the 
best-fit parameter values are Q.2,0 = 0-277^q[J^ anc [ n — 
0.70+| ;^. For the EoS of IDM, the best-fit value is w IDM = 
— 1.23^q?|. These results are consistent with the previous results 
of Basilakos & Plionis (2009) by using fewer observational H{z) 
data (n -0.30 and widm - -0.90, which favors the effective 
quintessence model with a constant EoS for which w > — 1). In 
Fig. [5] we show the constraint results from SNe la. By minimiz- 
ing the corresponding^ 2 , we find that the best-fit model values 




Fig. 5. The likelihood contours at the 68.3% and 95.4% confi- 
dence levels in the Q2,o - n and ^2,0 - widm planes provided by 
fitting model 2 to the SNe la data. 




Fig. 6. The likelihood contours at the 68.3% and 95.4% confi- 
dence levels in the Q2,o - n and ^2,0 - widm planes provided by 
fitting model 2 to the H(z)+SNe data. 

are Q 2 ,o = 0.288!°;}°* and n = 0.17!};}° (w lDM = -1.06^). 
For the combined data H(z)+ SNe shown in Fig. [6] we find 
the best-fit values are Q 2 , = 0.28 1!° °^ 2 and n = 0.22!° ^ 
(widm = -J --07!°)23)- Obviously, all of the above constraints are 
consistent with the concordance ACDM model and appear to be 
most consistent with the effective phantom model and a constant 
EoS for which w < — 1. 



5. Conclusions 

We have investigated the interacting dark matter (IDM) sce- 
nario mimicking either the ACDM model or the wCDM model, 
which can create the cosmic acceleration without dark energy 
(Basilakos & Plionis 2009). In our work, the scale of the effec- 
tive annihilation term kC\ and therefore the mass of DE parti- 
cles M x are constrained with different newly revised observa- 
tional data including H(z) and Union2 SNe la data. When mim- 
icking a A cosmology and using three different data combina- 
tions of H(z), SNe la, and SNe la+H(z), we have found that 
kC\ is quite small, which is consistent with the previous results 
in Basilakos & Plionis (2009), which used fewer observational 
H(z) data. Meanwhile, for the combined data sets, we obtain a 
more stringent upper limit to the effective annihilation term with 
log(A:Ci ■ Gyr) « -3.4. By relating the range of values of kC\ to 
the mass of the DM particle, we have inferred an apparent lower 
limit of M x ~ 10~ 8 6 Gev. Furthermore, when mimicking wCDM 
model and assuming that the particle creation term dominates 
(k = 0), we obtained the effective equation of state of IDM is 
consistent with the concordance ACDM model and appears to 
be consistent with the effective phantom model with a constant 
EoS for which w < — 1. 
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To sum up, we conclude that the interacting dark matter 
(IDM) model may provide a practical alternative to Dark Energy 
in explaining the present cosmic acceleration. The hope of prov- 
ing this conclusion should be pinned on analyses of future ob- 
servational data such as high redshift SNe la data from SNAP 
(Albrecht et al. 2006), more precise CMB data from the ESA 
Planck satellite (Balbi 2007), and complementary data, such as 
the X-ray gas mass fraction in clusters (Allen et al. 2004; Allen 
et al. 2008; Ettori et al. 2009), gravitational lensing data (Zhu 
1998; Sereno 2002), as well as gamma-ray bursts (GRBs) at 
high redshift (Schaefer 2007; Liang et al. 2008; Basilakos & 
Perivolaropoulos 2008; Liang & Zhang 2008; Wang & Liang 
2010; Gao et al. 2010; Liang et al. 201 1). 
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